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ABSTRACT 



HYDROFOIL SHIP 
CONSTANT LIFT CONTROL SYSTEM 
PART I: HYDRODYNAMICS 

MICHAEL ROY TERRY 
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NAVAL ARCHITECTURE and MARINE ENGINEERING 
and MECHANICAL ENGINEERING on 
23 May 1969, In partial 
fulfillment of the requirements 
for the Degrees of NAVAL ENGINEER, 
and MASTER of SCIENCE in MECHANICAL ENGINEERING. 



This paper analyzes the unsteady hydrodynamics of 
1) FULL INCIDENCE; 2) TRAILING EDGE FLAP; and 3) TAB FOIL 
hydrodynamic control devices. The device is to be used in a 
constant lift, wave alleviation control system on an OPEN- 
OCEAN HYDROFOIL SHIP. Each device is compared in POWER, DRAG, 
and WEIGHT. The conclusions reached are: 

1. Wave Alleviation is an important aspect of 
hydrofoil ship control system design. 

2. Unsteady Analysis is vital to understanding 
the response of a foil in waves. 

3. Negative Control System Power is required for 
wave alleviation by a tab foil with a lagging 
tab . 

4. Trailing Edge Fla p requires less overall 
system power than a Full Incidence Foil . 

5. Hydrofoil Ship Operation Policy will be greatly 
affected by unsteady hydrodynamics. 

6. Unsteady Foil Response is not a unique function 
of encounter frequency. 



Thesis Supervisor: Damon Ellis Cummings 

Title: Assistant Professor of Naval Architecture 
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SYMBOLS 



a » NACA camber loading parameter 

a « a hinge point in % of b 

A « Foil area 

b = Semichord 

c = Chord 

c a 3 hinge point in % of b 

C » Designates circulatory 

C 2 = Lift coefficient L/( (1/2 )pV 2 A) 

C l0 = Lift coefficient at ideal angle of attack 

C(k) = Theodorsen 2-D unsteady lift function 
d « Source position 

D = Drag 

e = Designates encounter 

g a Gravitation constant ft 2 /sec 

h a Keel to foilborne waterline height 

i “ /_! 

I a “ Main foil inertla/ft span 

Ijj = piap inertia/ft span 

J n (k) = Bessel Function of order n 
k a Reduced frequency 

L = Li ft /ft span 

L » Length 

m = Source strength/ft span 

M = Moment /ft span 

M = Mass of the ship 
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NC 3 Denotes non-circulatory 

Q * Flow rate 

s = Span 

t 3 Time 

t 3 Thickness 

T 3 Period 

U 3 Ship Speed 

V 3 Wave Phase Velocity 

w = Denotes wave 

x = Distance in horizontal direction 

Xo 3 Longitudinal distance from reference point to 

center of gravity 

y 3 Distance in vertical direction 

z 3 Distance traveled in vertical direction 

z 3 Vertical velocity 

•• 

z 3 Vertical acceleration 

Z 3 Vertical force 

a = Pitch angle of incidence about ba 

a p 3 Pitch angle of incidence magnitude 
a 3 Incidence angular velocity 

a 3 Incidence angular acceleration 

8 3 Pitch angle of flap about be 

8 0 3 Pitch angle of flap magnitude 

• 

8 3 Flap angular velocity 

8 “ Flap angular acceleration 
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A = Ship displacement 

A = Increment of change 

e = Small increment 

Hp = Propulsive coefficient 

n = Wave height 

ho “ Wave height magnitude 

n «* Wave orbital velocity 

fi » Wave orbital acceleration 

•• 

0 * Ship pitch acceleration 

X = Wave length 

it = 3 . 1 M 6... 

p » Fluid density in slugs/ft* 

4) = phase angle 

4 >i = lol 

*w s AH 

A4> = Phase angle between a, B 

= Angle of attack of ship to waves 
ip *» Stream function 

w « Circular frequency radians/sec 
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Renewed interest in the use of hydrofoil lifting sur- 
faces to raise a ship above the surface of the water 
occurred in the 1950's. Since that time, development of 
complete hydrofoil ship systems has progressed at a rapid 
rate. During the recent development stages, the encounter 
between ships and open-ocean waves has identified the pro- 
blem of ship motion control as being critical. The foil 
encounters lift changes caused by wave particle orbital 
motion. The wave particle orbital velocity is superposed 
on the uniform in-flow velocity (ship's speed through the 
water). The superposition results in an angle of attack 
change relative to the foil and, therefore, changes its 
lift. The changes in lift are disturbing forces to the 
ship and occur at the encounter frequency of the ship meeting 
the waves. Figure (1) depicts this situation. 



lb 



crest to trough 




This paper will address only open-ocean waves and ocean- 
going hydrofoil ships. 

During the early stages of development, attention was 
focused on the choice of foil configuration for best results. 
From that lively discussion and experimentation, two major 
design configurations developed. The first is commonly 
referred to as surface piercing , and the second as f ul ly 
submerged. Figure (2) depicts these configurations. 
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This paper will address only the fully submerged configura- 



tion . Hydrofoil ships with fully submerged foils are Inher- 
ently unstable in heave, pitch, and roll when the hull is 
out of the water, because the foil has very little sensiti- 
vity to depth below the water surface. Therefore, a control 
system must be provided. 

While foilborne, two modes of operation are possible. 
The first is the platform mode where the ship’s motion is 
maintained at constant height above the mean water level, 
at constant attitude about the height, and with no accelera- 
tions. The platform mode is only possible in wave heights 
that do not exceed the keel-to-foilborne waterline distance. 
The second is the contour mode where the ship’s motion is 
maintained at constant height above the local water surface 
at two points along the ship’s length. The contour mode 
is only necessary in wave heights that exceed the keel-to- 
foilborne waterline distance. Figure (3) depicts these 
modes . 





platforming 



figure 3 
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This paper will address only the platform mode and the wave 
systems that allow the use of that mode. 

The design of motion control for a hydrofoil ship must 
meet the three basic objectives outlined in Figure (^). 



STABILITY 



(Maintaining Commanded 
Equilibrium Attitude) 



MANEUVERABILITY 



(Changing to a New 
Commanded Equilibrium) 



WAVE ALLEVIATION 



(Reducing the Effect 
on Lift due to Waves) 



figure 4 

Most attention in the past has been concentrated on stability 
and maneuverability control functions. The effect of waves 
on the ship has been treated as a disturbance from equili- 
brium, and the stability elements of the control system have 
been 3ized to be able to correct the ship’s attitude and 
reduce the accelerations to a tolerable level. In contrast, 
this paper will address only the wave alleviation objective. 
The control system is usually composed of elements as 
shown in Figure (5). 
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SENSORS 



AUTOPILOT/COMPUTER 



SERVO 



HYDRODYNAMIC DEVICE 



(ELECTRONIC) 
(ELECTRONIC) 
(FLUID POWER) 
(LIFTING SURFACE) 



figure 5 

This paper will address only the Hydrodynamic Device* 
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II. PROJECT 



The objective of the above focus is now described. An 
open-ocean hydrofoil ship will operate, at times, when wave 
systems are present for the complete voyage. Operation in 
waves is, therefore, an important aspect of the design. Would 
the control system perform better than present day systems 
if the approach to its design considered only wave allevia- 
tion as the objective? The hypothesis 13 that the power 
level of the control system operating in waves could be 
significantly reduced if this approach was taken. Also, the 
wave alleviation control system could be integrated with a 
present day standard stability and maneuverability control 
system to produce a control system with overall performance 
superior to today's design. The stabilizing and maneuvering 
elements would effectively be operating as if they were in 
calm water, when in fact the hydrofoil ship is flying in 
open-ocean waves. Figure 6. 




figure 6 
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In the process of designing the optimum control system for 
wave alleviation, different solutions must be tested. The 
criteria best suited for this test is minimum overall system 
power. The three elements that contribute directly to over- 
all system power are shown in Figure (7)» 




Power to Operate 
the Control Surface 



Power to Overcome the Drag 
of the Hydrodynamic Device 



WEIGHT 



Power to Lift, Through the 
Lift/Drag Ratio, the Weight 
of the Control System 



figure 7 



For different designs, any one of the above elements of 
overall control system power may predominate. The task of 
designing the control system can be broken down into the 
design of the four basic elements — sensor, autopilot, servo, 
and hydrodynamic device as shown in Figure (5)» "Hydrofoil 
Ship Constant-Lift Control System: Part I Hydrodynamics" 

will address only the fourth element— Hydrodynamic Force 
Producing Device . 
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In summary. Figure (8) shows the steps in focusing this 



thesis within the overall project. 
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III. HYDRODYNAMIC DEVICE 



A. DEVICES 

The first step in the plan to find the minimum 
power device is to develop as many candidates as possi- 
ble for closer scrutiny. In general terms, the hydro- 
dynamic device must effect a force. Although there 
are a number of ways of effecting a force, this paper 
will consider only those that develop lift . This is a 
reasonable specialization because the hydrofoil ship 
utilizes lifting surfaces for lifting itself out of the 
water and for stability and maneuverability. 

Table (1) is an attempt at classifying the candi- 
dates . This paper will pursue only the following 
candidates : 

1. Full Incidence 

2. Trailing Edge Flap 

3. Tab Foil 

Appendix (A) describes the remaining candidates and 
discusses some of their characteristics in qualitative 
terms . 
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CAMBER AND ANGLE OP ATTACK CHANGERS 




Li. oil IP 



before proceeding with a quantitative analysis of the 
three candidates, a reference hydrofoil ship should be de- 
fined for Use as a common denominator for each design. 
Table (2) lists the parameters of the reference hydrofoil 
ship chosen. 



TABLE (2) 



Displacement 
Length Overall 
Design Speed 
Design Lift/Drag 



A 

L 

U 



300 Tons 
220 Ft 
40 kts 



L/D = 7 



Design Propulsive Coefficient n p " .6 
Each of 2 Main Lifting Foils NACA 16-408 
Thickness/Chord t/c = .08 

Mean Line a = 1.0 

Design Lift Coefficient C-^o “ .4 



Chord 
Semichord 
Span 
Weight 

Foil Inertia/Ft Span I 
1/4 Chord Trailing Edge Flap 
Chord 
Wt. 

Flap Inertia/Ft Span I 



c = 8.58 ft 

b - 4.29 ft 
s « 21.25 ft 
wt. = 7 tons 



= 140.0 Slug Ft 2 /ft 

= 2.145 ft 
a 1.75 tons 
« 8.80 Slug Ft 2 /ft 



Keel to Foilborne Waterline Height h = 14.0 ft 
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IV. POWER 



The first and nost important component of overall con- 
trol system power is the power to drive the hydrodynamic 
device through motions that will effect platforming in 
•waves. A reasonable model of motion in the vertical plane 
for the hydrofoil ship is 
• • 

Mz + Mx 0 6 * Z (waves) + Z(3hip + Z( Control (1) 

Motion) Motion) 

The effects of waves, ship motion, and control motion are 
modeled as linearly superimposed. For the platforming mode, 
the equation is reduced to 

Zero = Z( waves) + Z( control (2) 

motion ) 

The objective is to cancel the lift caused by waves by 
generating an equal and opposite lift caused by the control 
surface. From the required motion of the control surface 
and the resultant forces on it, the power required can be 
calculated 

Power » (motion velocity ) (summation of forces in ( 3 ) 

the direction of the velocity) 

The next task is to look at the effect of waves on a foil. 

A. QUASI-STEADY LIFT DUE TO WAVES 

From a qualitative, quasi-3teady point of view, the 
lift due to a wave passing a fully-submerged hydrofoil is 
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caused by the orbital velocity of the water particles super- 



posed with the velocity of the foil through the 'water. For 
a steady foil, the orbital velocity appears as an angle of 
attack change for the incoming flow. Figure (9). 



For the hydrofoil ship, the limiting wave condition for 
continued platforming is— total amplitude of the wave equal 
to the keel-to-foilborne waterline height. For the parti- 
cular hydrofoil ship chosen, this means a wave half-height 
of 7.C feet. For an open-ocean wave, the common wave length 
for this amplitude is 280.0 feet. At the ocean surface, at 
a given point, a reasonable model of the wave vertical posi- 
tion is 






figure 9 



Re{n 0 exp(-i2irV/A-t)} 



For a deep water x/ave 



V. =» /gX/2n 



27 



or 



n c Re{n 0 exp(-iV2irc/A • t)> 



( 6 ) 



The wave vertical particle velocity is then 

= Re{ri 0 (-i/2irg/X)exp(-i/2Trg/A • t)} (7) 

The amplitude of the orbital velocity is then 

| n | = no /2 tt g/l (8) 

For the chosen wave 

f n | “ 5.85 ft/sec (9) 

With the hydrofoil ship'3 velocity at 67.5 ft/sec, then the 
maximum angle of attack change due to vertical orbital velo- 
city is 

Tan a = 5.95/67.5 (10) 

or 

a « 5.06° (11) 

The maximum inflow velocity change due to horizontal orbital 
velocity is 



U x = 67.5 ± 5.95 ft/sec (12) 

To determine the effect on lift, we must look at the lift 
coefficient C^, For the 5.06° angle of attack change , flat 
plate model of the foil gives 



28 



( 13 ) 



ac^/aa = 2n 

or 

AC^ e 2nAa (14) 

AC X = .5^6 (15) 

so 

C lX « + .546 * .946 (16) 

For the hydrofoil ship chosen, the ratio of the lifts is then 

Lj /L 0 - 2.36 (17) 

In contrast, for* the 5.95 ft /sec Inflow velocity change , the 
ratio of the lifts is 



L t /L # = (Ui/U ,,) 2 



( 18 ) 



L,/L 0 « 1.15 



(19) 



From the above, the effect on lift due to the vertical compo- 
nent of orbital velocity is far greater* than that due to the 
horizontal component of orbital velocity. This paper will 
only consider the effect of the vertical component of orbital 
velocity. 

The use of quasi-steady analysis is not adequate for 
the case of the hydrofoil ship chosen. The encounter 
frequency, semichord, and ship's velocity combine to give 
a reduced frequency 



k =u e b/U 



( 20 ) 
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in the range where unsteady effects are most important, 
as will be shown in Section (B). Figure (10) describes 




u e = 2n/T e = 2nU e /L e (21) 



where 

L e » X/Cos^ (22) 

and 



U e = (V/Cos<|>) + u 



( 23 ) 



so 



u> e = (2ir/X)(V + UCosi}/) 



( 2 ^) 



but 

V 2 = gX/2Ti 



( 5 ) 



therefore 



X - 2irg/u) w 2 



( 26 ) 
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or 



s w u + 10 2 * UCosiJ//g (27) 

^ W Vil 

For the wave and ship chosen 

» e « .85 + .72(67. 5) (l)/32.2 = 2.36 Rad/Sec (28) 

Semichord b « A* . 29 ft 
Speed U = 67.5 ft/3ec 

so 



k = 2.36(4. 20 ) /Cl . 5 = .15 (29) 

or k is in the most critical range for unsteady effects as 
shown In Section B. 




figure 11 
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The following i3 an adaptation of Chapter 5 of 
Reference (1) for this particular problem. Please refer 
to Table ( 3 ) in this paper for the following discussion, 

A model of a wave passing a steady foil, shown in Figure (11), 
is expressed by Equation (30). The disturbing upwash velo- 
city UPW is equal to a complex phasor with magnitude — |n|, 
phase angle (relative to wave height) — <|> w , and being driven 
around the phasor diagram with angular velocity— oj q , 

Assuming negligible effects from free surface, negli- 
gible energy radiation on the free surface; idealized 2-D 
flat plate theory, with no separation or cavitation, the 
unsteady lift is expressed as the quasi-steady lift modified 
by the Theodorsen Function and additional Bessel Functions. 
These functions are expressed in Equations (31 to 3 * 0 . The 
Theodorsen Function is shown in figure (12). 




From Figure (12), it can be 3een that the Theodorsen Function 
reduces in magnitude (from 1 to 1/2) and has its maximum 
lagging phase angle change near k = .15; — exactly where the 
hydrofoil ship is operating. 
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CTkT = H^UkJ/Cll/^Ck) + ill/ 2 '(k)] (31) 
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The expression for the complex phasor LI Fl'W is shown 
in Equation (35) • The complex phasor moment MOMTW due to 
LiFTW acting about point ba (Figure (11)) is expressed in 
Equation (36). The complex phasor moment MOMTBW due to lift 
on the flap acting about point be (Figure (11)) is expressed 
in Equation (37)* These functions are best visualized with 
a phasor diagram. Figure (13) shows a phasor diagram scaled 
to illustrate the relationship between UPW, LIFTW, MOMTW, 
and MOMTBW for w e = ( + ) . 

NOTE: 

1. The diagram rotates clockwise . 

2. The uowa3h velocity leads the wave height by 90°. 

3. LIFTW lags UPW. 

4. MbMTW is in phase with LIFW. 

5. M0MT3W is l30° out of phase with LIFTW. 

6. LiFTW and MOMTW would have been larger in magnitude 
and in phase with UPv/ , if k = 0. 

Figure (l 1 *) shows a phasor diagram again scaled to illustrate 
the same relationships for <u e = (-). 

NOTE : 

1. The diagram rotates count ere lockwlse . 

2. The upwash velocity lags the wave height by 90°. 

3. Again LIFTW and MOMTW lag UPW. 

Jj. MOMTBW is 180° out of phase with LIFTV/. 

Figure (15) illustrates the meaning of u> e = (+,-). 
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astern sea 




1. Alpha Motion 

The following is again an adaptation of Chapter 5 of 
Reference (1). Please refer to Table (M) in this paper for 
the following discussion. A model of full incidence motion 
(Alpha Motion) (Figure ( 16 ) ) is expressed in Equations (38 
to MO) . It is expressed as a complex phasor with phase 
angle relative to the wave height— *<J> j , and being driven 
around the phasor diagram with angular velocity— w e . 

Again, assuming negligible effects from free surface, 
negligible energy radiation on the free surface; idealized 
2-D flat plate theory, with no separation or cavitation, 
the unsteady complex phasor lift (lTfTJ) is expressed in 
non-clrculatory velocity and acceleration terms, and circu- 
latory position and velocity terms modified by the Theodcrsen 
Function. YJWK is expressed in Equation (Ml). The unsteady 
complex phasor moment MOMTA due to XTIWa acting about point 
ba (Figure (16)) is expressed in Equation (M2). The unsteady 



3? 




38 



-2pUb 2 { ( 1/2) [Cos" 1 c - c/l-c k ] + [(l+c/2)/l-c 1 

- (c+l/2)Cos" 1 c]cTk7}*{U* ALPHA + b(l/2 - a) -ALPHADOT) C (43) 
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complex phasor moment MOMTBA due to lift on the flap acting 
about point be (Figure ( 16 ) ) is expressed in Equation (^3). 
Again, these functions are best visualized with a phasor 
diagram. Figure (17) shows a phasor diagram scaled to 
illustrate the relationship between AEMlA, AlMADbT, ALPHA- 
EdS 1 ?, LI FT A, MO.MTA, and MOMTBA for w e « ( + ). 

NOTE: 

1. The diagram rotates clockwise . 

2. ALPHADOT leads ALPHA by 90°. 

3. ALPHADDOT leads ALPHADOT by 90°. 

4. LI FT A lags ALPHA. 

9. MOMTA and MOMTBA lag LIFTA. 

6. LiFl’A will be in phase with ALPHA for k « 0 
because although LIFTA has contributions from 
ALPHADOT and ALPHADDOT, these are zero for k = 0. 
Figure (18) shows a phasor diagram to illustrate the same 
relationships with u) e = (-) . 

NOTE : 

1. Tne diagram rotates counterclockwise . 

2. ALPHADOT leads ALPHA by 90°. 

3. ALPHADDOT leads ALPHADOT by 90°. 

H. LIFTA lags ALPHA. 

2. Beta Motion 

The following is again an adaptation of Chapter 5 of 
Reference (1). Please refer to Table (5) in this paper for 



- T 




the following discussion. A model of trailing edge flap 
motion (Beta Motion) (Figure (. 16 )) is expressed in Equations 
(M to 1 l6). It is expressed in terms of the complex phasor 
ALFITa. The ratio of magnitude between BETA and ALPHA is — 
go/cto; the phase difference between BE^TA and ALPHA (+ for 
BEtfA leading) is — A$; and the phasor is being driven around 
the phasor diagram with angular velocity — w e . 

Again, assuming negligible effects from free surface, 
negligible energy radiation on the free surface, idealized 
2-D flat plate theory; with no separation or cavitation, 
the unsteady complex phasor lift (tlWB) is expressed in 
non-circulatory velocity and acceleration terms, and circula- 
tory position and velocity terms modified by the Theodorsen 
Function. LIFTB is expressed in Equation (47). The unsteady 
complex phasor moment MOM^B due to tilFTB acting about point 
ba (Figure (16)) is expressed in Equation (47). The unsteady 
complex phasor moment MbM^BB due to lift on the flap acting 
about point be (Figure (16)) is expressed in Equation (49). 
Again, these functions are best visualized with a phasor 
diagram. Figure (19) shows a phasor diagram scaled to 
illustrate the relationship between BETA, BETADOT, BET AD DOT, 
lTTO, ran, and ram for (o e - ( + ). 

NOTE: Same as Figure (17) • 

Figure (20) shows a phasor diagram to illustrate the same 
relationships with w e » (-). 

NOTE: Same as Figure (18). 
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D. REQUIRED FOIL MOTION 



Mode 1 «* ALPHA only (Full Incidence) 

Mode 2 = BETA only (Trailing Edge Flap) 

Mode 3 - ALPHA and BETA (Tab Foil) 

The wave Is given. The lift due to the wave can be 
solved explicitly by Section B. However, the required con- 
trol surface motion ( AtPilA and/or BCTS) is only implicitly 
expressed in Equations (41) and (47) of Section C. Both 
Equations (41) and (47) are complex expressions. Each one 
represents two equations either as the real and imaginary 
parts or the magnitude and phase angle. 

The strategy used to solve for TTEPTIa and/or BE^A is: 

1. Substitute in Equations (4l) and/or (47) expressions 
for ALPHA and/or BETA and their derivatives in terms 
of an unknown ct 0 and <|>|. (B 0 /a 0 and A<J> known con- 
stants ) 

2. Expand multiple angle expressions. 

3. Separate the real and imaginary parts of the modi- 
fied Equations (41) and/or (47). (-ReLlFTW, 

-ImLiFT'W) 

4. Divide out a 0 , Cos<J>j,p, and b. (-ReLlFTW/( aoCos^ pb) * 

5. Group the remaining terms into terms containing 
Tan 4>i— Xj, and constant terms — Y 2 , X 2 . 



-ReLlP i fW/(a e Cos4> l pb) * Tan4»i(Yi) + (Y 2 ) 



(50) 



-ImnTW/(aoCos<j)ipb) - Tan^Uj) + (X 2 ) (51) 

6. Divide the real equation by the imaginary equation. 

- ReUFW/ ( -I mYJWN) « R 

» (Tan<(» 1 (Y 1 )+(Y 2 ))/(Tan« 1 (X l ) + (X 2 )) (52) 

7. Separate out Tan 4>i 

Tan <j>, = (Y 2 -RX 2 )/(RX 1 -Y 1 ) (53) 

8. Solve for 4>i . 

9. Magnitude oo is 

a 0 » -ReLIP^W/[Cos4> 1 pb(Tan(J) 1 (Y 1 ) + (Y 2 ))] (5 1 *) 

10. The problem is solved. 

See Appendix C for details. With ALPHA and BETA, the equa- 
tions in Section C will give lift and moment from the motion. 

E. SUMMATION OP HYDRODYNAMIC FORCES 

Summation of LiFTW, LlFTA and/or L±FTB should equal 
zero. Summation of MbrWw , MOMTA and/or m6Mtb gives the 
total hydrodynamic moment about ba. Summation of MOMTBW, 
M< ^ MTBA and/or MbMTBB gives the total hydrodynamic moment 
acting on the flap about be. 



F. REQUIRED APPLIED MOMENTS 



The control system must apply a moment (MOMAP) on the 
control surface to produce the required motion. The dyna- 
mics of this interaction is expressed in Equations (55) 
and (56). 

I a *AEFHADDOf = MOMTW + MOMTA + MOMTB + MOMAPA (55) 

I b -BETADDOT = MoMtfBtt + tfOMfBA + MOMTBB + MOMAPB ( 56 ) 

From these equations, the required applied moments can be 
found. Figure (21) shows these relationships in phasor form. 




G. AVERAGE POWER 



The average power (PAVE) required to produce the 
motion is expressed in Equations (57) and (58). 

PAVEA => | ALPIIADOTl • |MOMAPA| ‘POWER FACTOR/2 (57) 

PAVEB = | WFKDQT | • | MOMAPB | • POWER FACTOR/2 (58) 

where the power factor is the cosine of the phase angle 
between the motion velocity and the applied moment. 

H. PROGRAM 

Two programs were written in Fortran IV for use at 
the Information Processing Center on the "IBM System 360" . 

The first program writes out the details of the power 
calculation, and the second writes out a minimum of identi- 
fying information and plots the power results on a "Calcomp 
Plotter”. An outline of the program is shown on Figure (22). 

The encounter frequency is calculated using Equation 
(27) as shown in Figure (10). The reduced frequency is 
calculated using Equation (20) but taking the absolute 
magnitude for use with Bessel Functions. The Theodorsen 
Function is calculated using Equation (31) by calling the 
"IBM Scientific Subroutine Package" routines BESJ and 
BESY . Because of the iterative method of calculation used 
in BESY, the reduced frequency parameter k can not be 
less than 0.05. This is guarded against in the program 
but results in errors in the Theodorsen Function for reduced 
frequency below 0.05. 





figure 22 
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The WAVLIF routine calculates the complex phasor 
forces generated by the wave using the analysis of Section 
B. The lift due to waves (LI i FTH*I ) i3 brought into the CONTRL 
routine. Here the complex phasor ALPHA and/or BETA are 
calculated to produce an equal and opposite lift to cancel 
ETFTw. The procedure used is described in Section D. The 
motion (ALPHA and/or BETA) i3 brought into the two routines 
LIFA and LIFB to calculate the lift caused by the motion, 1) 
as a check on the results of CONTRL; and 2) in order to 
calculate the moments about ba and be, using the equations 
of Section C. 

So far, the control surface motion and the hydrodynamic 
moments are known. The hydrodynamic moments due to waves 
and control surface motion are summed as in Section E. Next, 
the required applied moments are calculated as in Section F. 
The average power is calculated using the analysis of Section 
G. Note all results are 2-D per foot of span. An important 
requirement is the compatibility between mode specification 
and the input data values for BA = 0«/ao, and DPHE * A<j>. 

For Mode 1 ALPHA motion, 

BA « 0.0 DPHE ■ 0.0 
For Mode 2 PEl^A motion, 

BA * 1.0 DPHE » 0.0 
For Mode 3 AL^liA and tJE'TA motion 

DPHE / 0.0 
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The plotting routines plot the motion variable magni- 



tude and phase angle, and the average power as functions 
of encounter frequency. One plot per foil , the plotter 
treats each wave run as a point and connects each point 
with a straight line. The choice of wave input data should 
be arranged to give a somewhat continuous variation as was 
done for the series discussed in Section I. 

Appendix D contains a listing of the programs and an 
example of their written output. Examples of the plotted 
output can be found in Section I. 

I. RESULTS 

The two pieces of input data required for each run of 
the program are the foil data and the wave data. Three sets 
of calculations were performed for the three candidate foils, 
as described by Table (2): 

Mode 1 (Full Incidence) 

Mode 2 (Trailing Edge Flap) 

Mode 3 (Tab Foil) 

The first set is a series of waves with ship speed con- 
stant, wave length constant, and angle of attack (SI) between 
the ship and the waves varying from ’’dead ahead” to "dead 
astern”. Table (6) displays the wave parameters for each 
run of thi3 set. Figures (23), (24), and (25) are the compu- 
ter plotted results for the three foils run with this series. 
The series can be pictured as a slow 180° turn being executed 
by the ship at constant speed in regular waves of constant 
wave length. 
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RUN 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 



TABLE (6) 



SHIP ANGLE OF ATTACK TO THE WAVE, VARIATION 



HW“n 


WL-A 


u=u 


31=4/ 


PHEW=«}>w 


(FT) 


(FT) 


(Fr/SEC) 


(RADS) 


( RADS ) 


7.0 


280.0 


80.0 


0.0 


-1.5708 


7.0 


280.0 


80.0 


.3927 


-1.5708 


7.0 


280.0 


80.0 


.7854 


-1.5708 


7.0 


280.0 


80.0 


1.1781 


-1.5708 


7.0 


280.0 


80.0 


1.3745 


-1.5708 


7.0 


280.0 


80.0 


1.5708 


-1.5708 


7.0 


280.0 


80.0 


1.7672 


-1.5708 


7.0 


280.0 


80.0 


1.9635 


-1.5708 


7.0 


280.0 


80.0 


2.1599 


-1.5708 


7.0 


280.0 


80.0 


2.3562 


-1.5708 


7.0 


280.0 


80.0 


2.5526 


-1.5708 


7.0 


280.0 


80.0 


2.7489 


-1.5708 


7.0 


280.0 


80.0 


2.9453 


-1.5708 


7.0 


280.0 


80.0 


3.1416 


-1.5708 



52 



The second set is a series of waves with ship speed (U) 



varying from 90 ft /sec down to 25 ft/sec, wave length con- 
stant, and angle of attack between the ship and the waves at 
"dead ahead" and "dead astern". Table (7) displays the wave 
parameters for each run of this set. Figures (26), (27), 
and (28) are the computer plotted results for the three foils 
run with this series. The series can be pictured as a slow 
deceleration by the ship in ahead regular waves of constant 
wave length and a slow acceleration by the ship in astern 
regular waves of constant wave length. 

The third set is a series of wave 3 with constant ship 
speed, wave length (WL) varying from 280 ft to 140 ft, and 
angle of attack between the ship and the waves at "dead 
ahead" and "dead astern". Table (8) displays the wave para- 
meters for each run of this set. Figures (29)* (30), and 
(31) are the computer plotted results for the three foils 
run with this series. The series can be pictured as a slow 
transit by the ship at constant speed in ahead regular waves 
of reducing wave length and a slow transit in astern regular 
waves of increasing wave length. 

Additional trial and error runs were made in order to 
discover a useful motion ratio 0 o /a o (BA) and a useful phase 
angle difference A<J> ( DPHE) for Mode 3 (Tab Foil). The search 
is not complete, but BA » 2.0 and DPIIE ® -.3927 radians 
appeared best and were utilized for the Mode 3 runs. 
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RUN 

1 

2 

3 

It 

5 

6 

7 

8 

9 

10 

11 

12 

13 

l^t 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 



TABLE ( 7 ) 



HW=n 

( FT ) 

7.0 

7.0 

7.0 

7.0 

7.0 

7.0 

7.0 

7.0 

7.0 

7.0 

7.0 

7.0 

7.0 

7.0 

7.0 

7.0 

7.0 

7.0 

7.0 

7.0 

7.0 

7.0 

7.0 

7.0 

7.0 

7.0 

7.0 

7.0 



SPEED VARIATION 



WL=A 


U«U 


( FT ) 


( FT/SE 


280.0 


90.0 


280.0 


85.0 


280.0 


80.0 


280.0 


75.0 


280.0 


70.0 


280.0 


65.0 


280.0 


60.0 


280.0 


55.0 


280.0 


50.0 


280.0 


45.0 


280.0 


40.0 


280.0 


35.0 


280.0 


30.0 


280.0 


25.0 


280.0 


25.0 


280.0 


30.0 


280.0 


35.0 


280.0 


40.0 


280.0 


45.0 


280.0 


50.0 


280.0 


55.0 


280.0 


60.0 


280.0 


65.0 


280.0 


70.0 


280.0 


75.0 


280.0 


80.0 


280.0 


85.0 


280.0 


90.0 



SI “\|) 


PHEW =< f>w 


( RADS ) 


( RADS ) 


0.0 


- 1.5708 


0.0 


- 1.5708 


0.0 


- 1.5708 


0.0 


- 1.5708 


0.0 


- 1.5708 


0.0 


- 1.5708 


0.0 


- 1.5708 


0.0 


- 1.5708 


0.0 


- 1.5708 


0.0 


- 1.5708 


0.0 


- 1.5708 


0.0 


- 1.5708 


0.0 


- 1.5708 


0.0 


- 1.5708 


3.1416 


- 1.5708 


3.1416 


- 1.5708 


3.1416 


- 1.5708 


3.1416 


- 1.5708 


3.1416 


- 1.5708 


3.1416 


- 1.5708 


3.1416 


- 1.5708 


3.1416 


- 1.5708 


3.1416 


- 1.5708 


3.1416 


- 1.5708 


3.1416 


- 1.5708 


3.1416 


- 1.5708 


3.1416 


- 1.5708 


3.1416 


- 1.5708 
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TABLE (8) 



WAVE LENGTH VARIATION 





HW«n 


WL»A 


U=U 


SI*=i}> 


PHEW=4>w 


RUN # 


(FT) 


(FT) 


(FT/SEC) 


(RADS) 


( RADS ) 


1 


3.5 


140.0 


80.0 


0.0 


-1.5708 


2 


4.2 


168.0 


80.0 


0.0 


-1.5708 


3 


4.9 


196.0 


80.0 


0.0 


-1.5708 


4 


5.6 


224.0 


80.0 


0.0 


-1.5708 


5 


6.3 


252.0 


80.0 


0.0 


-1.5708 


6 


7.0 


280.0 


80.0 


0.0 


-1.5708 


7 


7.0 


280.0 


80.0 


3.1416 


-1.5708 


8 


6.3 


252.0 


80.0 


3.1416 


-1.5708 


9 


5.6 


224.0 


80.0 


3.1416 


-1.5708 


10 


4.9 


196.0 


80.0 


3.1416 


-1.5708 


11 


4.2 


168.0 


80.0 


3.1416 


-1.5708 


12 


3.5 


140.0 


80.0 


3.1416 


-1.5708 
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figure 23 
Full 
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figure 24 
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figure 27 
Trailing 
Edge Flap 
(U) Var. 

B = l/3l 
Y - LP 
Q = paveb 
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figure 28 

Tab 

Foil 

(U) Var 
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figure 29 
Full 

Incidence 
(WL) Van 

A = I ai 
Z = la 
P = pavea 
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figure 30 
Trailing 
Edge Flap 
(WL) Var. 
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figure 31 

Tab 

Foil 

(WL) Var. 

A = I a I 
\P\ = 2 • la, | 

Z =/a 

Lp = /a -22 '4 
P - pavea 
Q = - paveb 



All the results are in error for the range of reduced 
frequency, -.05<k<+.05 , because of the flx-up procedure 
mentioned In Section H. The error Is dramatic because of 
the sensitivity of power to the related phase angles of its 
components. The error is zero at k = 0, and beyond the 
range k =±.05. The range for k translates Into a differ- 
ent range of u> e for each of the three wave series. The 

reason is that k is composed of both u> and U, which are 

6 

not solely functions of one another. Equation (27). For 
a given U, co e remains dependent on X and The following 
are the ranges of w e which are in error for: 

Series 1. (SI) Variation 

-.85<w <+.85 (Rads/Sec) 

Series 2. (U) Variation 

-.94<a) e <+.29 (Rads/Sec) 

Series 3. (WL) Variation 

Runs not In the range 

The plotted results are somewhat misleading in a part 
of the u> e range, because the plotter plots straight lines 
between points. The following are the ranges In u> e where 
there are significant gaps in plotted points: 

1. (SI) Variation 



None 



2. (U) Variation 

+ .28<u> e <+l. Ml (Rads/Sec) 

3. (WL) Variation 

-.9M«d <+2.6M (Rads/Sec) 

e 

Aside from the fix-up error and taking note of the range for 
the gap in points, the plotted results give a reasonably 
good picture of the power required. 

In the case of Mode 3 ALPHA and BETA Motion (Tab Foil), 
the plotter plots the Beta power as Q a - PAVED. This opera- 
tion places the Q result relative to the P result (PAVEA) 
such that the distance between them is the total power re- 
quired by the system. If the distance between Q + P is in 
the positive direction, the total power required is positive . 
If the distance between Q P is in the negative direction, 
the total power required is negative . The above is true 
regardless of where the plots lie on the absolute scale. 

THE FOLLOWING ARE COMMENTS ON EACH FOIL RESULT: 

1. Alpha Motion (Full Incidence) Mode 1 

a. (SI) Variation* Figure (23) shows a decreasing 

average power requirement from "dead ahead" around to 
quartering seas where the encounter frequency goes to zero 
(power goes to zero); and then continuing around to "dead 
astern", the power increases again for w e = (-). The maximum 
power occurs at SI = 0 and it is positive. 
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The magnitude of the Alpha motion is constant, because 
the disturbance is unchanged in magnitude and the ship's 
(lift generating) speed Is constant. Note, at w * 0 
(SI » 112°), the magnitude of Alpha does not go to zero. 

This is because, depending upon what the phase angle of the 
wave is that the ship is running along (crest, back- face, 
trough, or front-face) , Alpha will have to be a constant 0, 
(+), 0 , (-). 

The phase angle of Alpha relative to the wave height 
changes from lagging 90° + e for u * ( + ), to leading 90 ° - e 
for w e » (-). At w e * 0, the phase angle is 90° from the 
wave height; but, more important, it is 180° out of phase 
with the wave upwash velocity UW as is expected from quasi- 
steady analysis. The basic characteristic of shifting from 
lagging to leading and being out of phase with UPW at u e s 0 
is common to all the foils In all the wave series. 

b. (U) Variation. Figure (26) shows increasing 

magnitude for the required Alpha motion, as speed drops from 
90 ft/sec in ahead seas (lesser w *» ( + )). The magnitude 
decreases as the ship speed Increases in astern seas (greater 
U) e B (-)). 

The phase angle goes through a similar variation as in 
the SI series. 

The average power decreases as the ship slows in ahead 
seas. In astern seas, the power dips down through zero at 
d) e = 0 and then climbs back up as speed increases. The gap 
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between points at ± 25 ft/sec is necessary because of the 
limited capability of the lifting foils to support the 
ship at low speed. This is evident by the rapidly increasing 
motion magnitude (shown in Figure (26)) required to produce 
the lift at low speeds. 

c. (WL) Variation. Figure (29) shows the magnitude 
of Alpha motion increases slightly as the wave length 
increases in ahead seas. This is because the increased 
wave length presents an increasing disturbance. The same 
variation appears in astern seas. 

The phase angle changes in a similar fashion as in the 
SI series. 

The average power decreases as wave length increases 
in ahead seas (lesser u> e *= ( + )). In astern seas, the power 
again decreases as wave length increases (lesser w e = (-)). 
Note, there is no effect of C(k) error in this plot, because 
<*> e is out of the range, 

2. Beta Motion (Trailing Edge Flap) Mode 2 

a. (SI) Variation, Figure (24) shows the magnitude 

of the Beta motion constant as in the Alpha case. The magni- 
tude of Beta is larger than the Alpha case which is to be 
expected for a 1/4 chord flap compared to full incidence. 

The phase angle goes through a similar variation as in 
the Alpha case. 
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b. (U) Variation. Figure (27) shows the magnitude 

increasing in a similar manner as in the Alpha case. Again, 
the magnitude of Beta is larger than that for Alpha. 

The phase angle change is similar to the (SI) variation. 

The average power follows a similar trend as in the 
Alpha case; but, again, it's sensitivity to error in C(k) , 
for the range -.9^<w e <+. 29, is much greater. The only point 
that is correct in this range is at w e = 0. Again, it is 
speculated that, if correct values for c(k) were used, the 
power curve would fair into zero at u> e = 0 and climb back 
up again in a similar fashion to the Alpha case. 

c. (WL) Variation. Figure (30) shows the magnitude 

of Beta increasing in a similar fashion as in the Alpha case. 

The phase angle changes in the same manner as in the 
Alpha case. 

The average power Increases for increasing wave length 
in ahead seas (lesser u e « (+)). This result is in contrast 
with the Alpha result where power decreased for increasing 
wave length in ahead seas. The average power in astern seas 
decreases for increasing wave length (lesser u> e ** (-)). 

Note, there is no effect of C(k) error in this plot, because 
u> e is out of the range. 

3. Alpha and Beta Motion (Tab Foil) Mode 3 

a. (SI) Variation. Figure (25) shows the magnitude 

of Alpha as a constant straight line as In the Alpha only 



73 



case. The associated Beta magnitude is two times the Alpha 
magnitude; however, it is not plotted. This is true for the 
(U) variation and (WL) variation as well. 

The Alpha phase angle has a similar trend, as before, 
but the w e » (±) variations are displaced away from 90°. 

The Alpha motion now lags the wave height by 90° - e for 
w e “ ( + )> and leads by 90° + e for u> e «* (-). The ramp-type 
Jog near u> e * 0 is caused by the plotter moving from w e = 

+ .16 *► -.14. The jog would approach a step type for points 
approaching w e «* 0. The Beta phase angle lags the Alpha 
phase angle by 22 1/2° for both u> e = (±); however, it is not 
plotted. This is true for the (U) variation and (WL) varia- 
tion as well. 

The results for power are important. As was pointed out 
earlier in this section; if the distance between Q (Beta 
power) -*■ P (Alpha power) is in the negative direction, the 
total power to drive the motion is negative (the control 
system must absorb power from the motion). Figure (25)'s 
plot of Q, P is exactly this case. Figure (32) demonstrates 
in qualitative form the relative position of the motion 
phasors and their lift phasors. 
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Xw. 




Figures ( 33 ) and ( 3*0 demonstrate in qualitative form how 
these phasors contribute to negative power. 
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b. (U) Variation. Figure (28) shows the magnitude 

of Alpha increasing as in the Alpha only case. 

The Alpha phase angle variation is similar to the 
Mode 3 SI case. 

The average total power is negative over the range of 
ui e except near w e * +1.50. Here the power required by Alpha 
is positive, but the power required by Beta is Ju3t enough 
negative for the total power to be zero. 

c. (WL) Variation. Figure (31) shows a slightly 

Increasing Alpha magnitude as in the Alpha only case. 

The Alpha phase angle is changing in a manner similar to 
the SI case. 
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The total average power starts out positive for small 



wave lengths in ahead seas. At u e = +4.28, the total power 
is zero and then goes negative for decreasing w e . The 
average negative power for Beta is not enough negative to 
cancel the positive Alpha power for w e +4.28. Note, the 
C(k) error does not effect this plot. 

From the above results, it is clear that an under- 
standing of the unsteady effects on foils is not only 
necessary but also rewarding. Table (9) li3ts some of the 
important power results. 



TABLE (9) 



(Full Incidence) 
Alpha Motion 
Maximum ± 
Average Power 
ft- lb/sec 
per ft span 



(Trailing Edge Flap) 
Beta Motion 
Maximum ± 
Average Power 
ft-lb/sec 
per ft span 



(Tab Foil) 

Alpha and Beta Motion 
Maximum ± 
Average Power 
ft-lb/sec 

per ft span 



SI 



+ 738.5 



+ 468.0 



- 154.9 



U 



+ 774.5 



+ 534.7 



- 153.5 



WL 



+ 1105.2 



+491.2 



+ 23.6, - 133.8 



V. DRAG 



The components of drag; for a foil are: 

Skin Friction 
Pressure Drag 
Induced Drag 
Wave Drag 
Interference Drag 

Each of these has a contribution from: 

Steady Motion 
Waves 

Foil Motion 

The task of sorting out all of these parts is formidable. 

For a 2-D analysis, the induced parts can be eliminated. 

For a small strut, the interference drag can be neglected. 

For a well-submerged foil, the wave drag can be neglected. 
Assuming small motions with no separation, cavitation or 
ventilation, the pressure drag can be neglected. 

This leaves the effects of steady motion, waves, and 
foil motion on skin- friction drag. The temptation is to 
linearly superpose these three effects as was done for lift. 
Yet the model for skin friction is the boundary layer. 

Reference (2) describes the effect on foil-skin friction 
by a uniform inflow with a harmonic perturbation superim- 
posed. This is a model of the effect of waves on the skin 
friction of a steady foil. The shearing stress lead3 the 
velocity perturbation (horizontal orbital velocity for the 
wave), and the phase angle approaches ^5° for k -*■ 00 . 
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An oscillating cylinder in a fluid at rest is discussed, 
and an oscillating infinite flat plate in a fluid at rest is 
discussed. The problem of a finite chord oscillating flat- 
plate foil in uniform flow is not discussed. No mention is 
made of the advisability of superposing such a solution on 
the effects due to waves. 

Although the required motion of the control surface 
results in zero total lift, the lift distribution 1 b differ- 
ent for the wave and for the motion. This is evidenced by 
the different values for wave moments and motion moments. 
Therefore, the pressure distribution will oscillate back 
and forth across the foil. It is speculated that the problem 
can no longer be solved by linear superposition, even if 
there was a solution to the foil motion effects. The solu- 
tion of the coupled effects of waves and motion appears 
beyond analytical solution today. 

All that is left is skin-friction drag from steady 
motion in uniform flow. The result for the three candidates, 
with the same overall foil geometry, is that there is no 
significant difference in skin-friction drag . 

This statement is adequate for comparing the overall 
power of the three candidates. However, if comparisons of 
these foils were being made with other systems, such as 
spoilers, the above rationalization is not adequate because 
ventilation and separation play a much different part in the 
drag of such systems. 



VI. WEIGHT 



Any statement about weight without a first iteration 
of the complete control system design must be qualitative. 
However, for the effect of weight to influence the choice 
of control system, an absolute figure must be found. For 
the purpose of emphasizing that power must be supplied to 
lift the weight of the control system, an approximate 
analysis will be pursued. 



A. REFERENCE FULL- INCIDENCE SYSTEM 

To use as a reference system, the weight of a full- 
incidence (Alpha motion) control system is estimated as 
follows : 



Group I Foil Weight 

Group II Mechanical Drive Links 

Hydraulic Servo Cylinders 
Hydraulic Fluids and Pipes 
Hydraulic Pumps and Drive 

Group III Autopilot and Computer 

Group IV Sensors 



=15,700 lbs 

= 650 lbs 

- 650 lbs 

» 1,070 lbs 

« *100 lb3 

■ 200 lbs 

* 200 lbs 



7 *io 


lbs/ft 


30 


lbs/ft 


30 


lbs/ft 


50 


lbs/ft 


20 


lbs/ft 


10 


lbs/ft 


10 


lbs/ft 



B. COMPARATIVE WEIGHT 

The Group I, III, and IV weights for the three candidates 
will not be significantly different. The Group II weights 
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will approximately vary linearly with average power. For 
comparison, from Table (9): 



Trailing Edge Flap Max Power 
Full Incidence Max Power 



53*,.Z 

1105.2 



.48 



Tab Foil Max Power « -154 .9 _ •, j. 

ji’ull Incidence Max Power 110$. 2 

(Here the maximum negative power is taken to size the 
servo machinery even though it must absorb the power.) 



C. TRAILING EDGE FLAP SYSTEM WEIGHT 

Group II » 70 lb 3/ft 

D. TAB FOIL SYSTEM WEIGHT 

Group II « 20 lbs/ft 



E. POWER 

The weight becomes power through the lift/drag ratio 
and the speed of the hydrofoil ship: 

Power » (Drag/Lift) ‘Weight ‘Speed (59) 

For comparative purposes, the lift/drag ratio is 7 and the 
speed is 80 ft/sec. Power results due to weight are as 
follows : 

Full Incidence =< 1500 ft-lb/sec per ft span 

Trailing Edge Flap * 800 ft-lb/sec per ft span 

Tab Foil e 230 ft-lb/sec per ft span 
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VII. CONCLUSIONS AND RECOMMENDATIONS 

In summary, the contributions of power, power due to 
drag, and power due to weight have been developed for three 
candidate hydrodynamic devices for use in a constant lift, 
wave alleviation, control system for a hydrofoil ship. 

A. OVERALL CONTROL SYSTEM POV/ER 

A servo system efficiency is required to convert the 
control surface power into input power to the system. An 
efficiency of .6 is assumed and has been applied to the 
maximum power results in Table (9)« 



Table 


(10) shows the results for overall 


control system 


power. 


TABLE 


(10) 






Pull Incidence 
H-ib/sec 
per ft span 


Trailing 
Edge Flap 
7t-lb7sec 
per ft span 


Tab Foil 
ft-lb/sec 
per ft span 


Power 


+ 1840. 


+ 890. 


+ 40. 


Drag 

Power 


— 


— 


— 


Weight 

Power 


+ 1500 . 


+ 800. 


+ 230. 


Total 


+ 3340. 


+ 1690. 


+ 270 . 



87 



From this analysis, the tab foil looks very promising for 
use as the hydrodynamic device in a constant lift control 
system. Development of the sensor, servo, and autopilot 
system designs for use with the tab foil is recommended. 

If the choice is restricted to either full incidence 
or trailing edge flap, the trailing edge flap looks most 
promising. The development of a control system utilizing 
a full incidence foil for stability and maneuverability , 
and a trailing edge flap for wave alleviation , is recommend- 
ed. 

But the most Important recommendations ares 

1. Wave alleviation should be an important considera- 
tion in the design of an open-ocean hydrofoil ship. 

2. Unsteady lift analysis should be used during the 
design of the control system for an open-ocean 
hydrofoil ship. 

B. SENSORS 

The control system must be able to sense the magnitude 
of the upwash velocity, the direction, and the encounter 
frequency of the wave the foil is about to see, enough in 
advance, in order to control proper foil motion. Wave height 
Information is not enough. It is suggested that wave height 
and orbital velocity be measured ahead and behind the ship. 

C. MOTION 

Figure (35) 3hows full incidence, trailing edge flap, 
and tab foil motion. The amplitudes have been exaggerated. 
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figure 35 



but their relative magnitudes have been maintained. The 



most interesting motion is that of the tab foil. With 
22 1/2° lagging phase angle for Beta (-A$), the Alpha motion 
is caused by the moments generated by the Beta motion. The 
motion appears "natural" in contrast to the 22 1/2° leading 
system shown below it. Figure ( 36 ) shows motions for 90° 
lagging, l80° out of phase, and 90 ° leading. 



£*- cOe "6 ^ 




Of these, only the 90° lagging appears "natural" • Preliminary 
checks with the computer program for 80 ft/sec speed in 
limiting ahead seas, indicate: 

90° Lagging } Alpha Power << 0 
Beta Power << 0 



90 



180° Out of Phase ) 



0 << Alpha Power 



Beta Power < 0 

90° Leading } 0 << Alpha Power 

0 << Beta Power 

The "natural” looking motions require negative power. The 
other motions require positive power. References (3) and 
(4) discuss the use of a tab to initiate aileron motion. The 
tab is moved in the opposite direction from the direction the 
aileron will move. At first, this sounds like the aileron 
tab should be 180° out of phase. However, the Beta motion 
must lead (in time) to cause the motion. This lead (in time) 
in the opposite direction is a lag relative to the aileron 
motion. 

For the hydrofoil ship, lags near 90° result in large 
negative powers to be absorbed (like the flutter situation). 
If small lags around 22 1/2° are used, small negative powers 
result. 

D. EFFECTS ON HYDROFOIL SHIP OPERATION 

In the limiting sea (ahead or astern), slowing down the 
ship will reduce the control system power but will increase 
the magnitude of the control surface motion. Since most 
foil arrangements have a mechanical limit to foil travel, 
slowing down may not be the best operational change. Running 
with the wave requires no foil motion and no control system 
power. Running with the crest or trough requires no control 
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surface angle. Running with the face or back requires full 
- or + control surface angle. Running with the crest is 
recommended because no control surface angle is required and 
maximum submergence of the foils is obtained. 

The angle of attack to the sea will reduce the control 
system power but not the foil motion magnitude. When running 
along the wave, the face or back of the wave present no span- 
wise flow; but the control surface will be at full - or + 
angle. Running along the crest or trough will present a 
small spanwise flow to the foils and struts, but no control 
surface angle will be required. Running along the trough is 
recommended because no control surface angle is required and 
maximum submergence of the foils, on the beam, is obtained. 

Finding smaller waves decreases control surface motion 
in all cases, and will reduce the control system power for 
a trailing edge flap system in ahead seas. But, smaller 
waves increase control system power for the full incidence 
system in both ahead and astern seas. Smaller waves Increase 
control system power for the trailing edge flap system in 
astern seas. 

It is recommended that wave alleviation unsteady effects 
be analyzed before determining operational policy. 

E. STATISTICAL ANALYSIS 

The foil responds to wave upwash velocity as an input. 
Therefore, an orbital velocity spectrum with direction infor- 
mation is required. This is a 3-dimensional-surface from 
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statistical wave data. Correspondingly, the foil response 
spectrum is no longer a single curve but is a 3-dimensional- 
surface for a given ship’s speed . The complete picture is a 
series of response surfaces for a range of speeds. 

It is recommended that an attempt be made at bringing 
the unsteady foil response analysis into a statistical for- 
mat. Prom this, a complete and realistic description of 
hydrofoil ship energy requirements could be made. 
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APPENDIX A: HYDRODYNAMIC DEVICES FOR SECTION III-A 

I. ANGLE OF ATTACK CHANGERS 

A. TRANSLATING FOIL (SUSPENDED FOIL) 

The foil vertical velocity must match that of the incom- 
ing wave. If the speed of the foil is zero, * w w . The 
vertical velocity of the incoming wave will be |n| *= n&> w . 

The velocity of the foil must be |z| » |n|; but, if u> e ® u> w , 
than |z| « | n | . In other words, the foil motion would exact- 
ly follow the wave surface. If speed is increased, the 
magnitude of w e Increases and |z| < |n|. Such large excur- 
sions would probably not be feasible. A patent for a similar 
scheme is in Reference (5)« ? 

II. CAMBER CHANGERS 

A. PARTIAL CHORD FLUID (JET FLAP) 

Much work has been done on jet flaps for high \ift on 
aircraft. The feasibility of pumping water through the foil 
has not been determined as yet. Fluidic controls could be 
employed to change the direction of the Jet forfwave allevia- 
tion. See References (6) and (7). 

B. UPPER SURFACE MECHANICAL (MEMBRANE) 

This scheme would employ a second skin for the .tipper 
surface of the foil. The space between could be pressurized 
to change the membrane shape. The upper surface lift load 
would have to be taken in tension in the membrane. The 
structure appears to be the greatest problem. 
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C. UPPER SURFACE FLUID (FLUID CAMBER) 

The concept of using a layer of fluid on the upper 
surface to change the apparent thickness to the oncoming 
flow was pursued to some detail early in the project. The 
following is an approximate analysis for such a scheme. 

In order to change the lift coefficient of the foil 
from .4 to .9, the camber must change from 2 , 2 % of chord to 
5 , 2 % of chord, or a change of 3$ of chord. In order to 
change the camber by 3*» the thickness must be changed by 6 % 
of chord. 

The proposed scheme of using "fluid camber" to change a 
foil's lift characteristic calls for pumping fluid from the 
leading edge on one surface of the foil and retrieving the 
same fluid back into the foil at the trailing edge. 
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After some fruitless searching for a potential flow 
model for the concentrated Jet at the leading edge, a half- 
plane source and sink were chosen as a simple model of both 
entry and exit of the flow. The stream function for a source 
and sink of equal strength m and separated by a distance d is: 



m 



rn — ** 1 -2y(l/2)d 

Tan xJVyJ-rr/2.d) a 



( 1 ) 
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The stream function for a uniform flow in the x direction 



at speed Uo is: 



* * U„y (2) 

For the superposition of these flows: 

* = U.y + fjTan" ( 3) 

We need the half plane thickness of the ellipse formed at 
■ 0 to be 6 % of chord, or y/d * .06. With “ 0, x ■ 0, 
y » .06d, we can solve for m. 



m » -(.06d)2v 

TJ7 ^an-TTTan H)' 



O) 



where 



= -2 ( . 06) ( l/2)d 

Tan(B) r;o6dTi-m72ldF 



Tan(B) = .2^3 



B = -166. 67° 



(5) 

( 6 ) 
(7) 



m = .1295* d*U 0 



( 8 ) 



But we need only 1/2 this flow but across the span of the 
foil. 



* .0648* Uo.ds 



(9) 



Q « .06^8*U 0 »A ft*/3ec 



( 10 ) 
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For an ocean hydrofoil ship of 300 tons at 40 knots 



with lift coefficient of .4, requires this area: 



A « 



L 



(ID 



(i/ , 2)pV2C 1 



A “ 




( 12 ) 



This gives us a required flow rate for the "fluid camber" 
system of: 



The "fluid camber" scheme is not feasible for the full scale 
hydrofoil ship. 

III. CAMBER AND ANGLE OF ATTACK CHANGERS 
A. ROTATING FOIL/PARTIAL CHORD FLUID (JET TAB) 

This scheme would utilize the Jet momentum and effect on 
the overall flow to force the Alpha motion. Reference (7) 
discusses such a scheme. Hopefully, the flow rate required 
in the Jet could be reduced from that in a Jet flap system. 

3. TRANSLATING FOIL/PARTIAL CHORD MECHANICAL (SUSPENDED FLAP) 
C. TRANSLATING FOIL/PARTIAL CHORD FLUID (SUSPENDED JET FLAP) 
Both schemes would use relative motion of a part of the 
chord to reduce the excursion magnitude of the foil. These 
schemes 3eem worth further analysis. 



Q » .0648 (67.5)068) « 1610 ft*/sec 

or = 7.22 x 10 s gal/min (13) 
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IV. LIFT SPOILERS 



A. MECHANICAL (SPOILERS) 

Mechanical Spoilers have been used in aircraft as 
ailerons. The scheme for a hydrofoil ship wave alleviation 
would require the foil to be at the angle of attack necessary 
to produce the maximum lift needed. Spoilers would then 
reduce that lift down to the minimum lift needed. The power 
and weight for such a system will be small, but the drag is 
expected to be great. 

B. FLUID (VENTILATION) 

The scheme of ventilating portions of the lift distribu- 
tion has been developed by SUPRAMAR. However, the same pro- 
blems of drag enter as in Section A. It may be that overall 
system power for such a scheme would be small compared to 
full incidence. 
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APPENDIX B 



DERIVATION OF MotWbW 
FOR SECTION IV-B 
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APPENDIX C 

DETAILS OP THE IMPLICIT SOLUTION 
FOR At?’:! A AND EeTTA FOR SECTION IV-D 
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APPENDIX D 

PROGRAM LISTINGS FOR SECTION IV-H 
I. COMPUTER PROGRAM FOILDYN (WRITING) 
II. COMPUTER PROGRAM FOILDYN (PLOTTING) 
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AVERAGE POWER REQUIRED FOR BETA® -35.54822 DUE TO POWER FACTOR® -0.02524 



